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Abstract 

We report an updated study of the suppressed decays B~ — > [K + ir~]£)K~ and B~ — > 
[K + ir~]£,n~ where [K + ir~]r, indicates that the K + ir~ pair originates from a neutral D meson. 
A data sample containing 386 million BB pairs recorded at the T(45) resonance with the Belle 
detector at the KEKB asymmetric e + e~ storage ring is used. This decay mode is sensitive to the 
CKM angle fo. We do not see a significant signal for B~ — > [K + it~]e)K~ , and we set a limit 
on the ratio of B decay amplitudes r# < 0.18 at the 90% confidence level. We measure the CP 
asymmetry of the B~ -> [K + TT-] D ir- mode, A D7r = 0.10 ± 0.22 (stat) ± 0.02 (syst). 
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INTRODUCTION 



Precise measurements of the elements of the Cabibbo-Kobayashi-Maskawa matrix [1] 
constrain the Standard Model and may reveal new physics. However, the extraction of 
the Unitarity Triangle angle 03 [2] is a challenging measurement even with modern high 
luminosity B factories. Several methods for measuring 3 use the interference between 
B~ — ► D°K~ and B~ — > D°K~, which occurs when D° and D° decay to common final 
states [3, 4]. CP violation occurs when both weak and strong phase differences between the 
amplitudes are non-trivial. As noted by Atwood, Dunietz and Soni (ADS) [5], CP violation 
effects are enhanced if the final state is chosen so that the interfering amplitudes have 
comparable magnitudes; the archetype uses B~ — > [K + ty~]dK~ , where [K + tc~]e, indicates 
that the K + tt~ pair originates from a neutral D meson. In this case, the colour-allowed 
B decay followed by the doubly Cabibbo-suppressed D decay interferes with the colour- 
suppressed B decay followed by the Cabibbo-allowed D decay (Fig. 1). Previous studies 
of this decay mode by BaBar [6] and Belle [7] have not found any significant signals for 
B~ — y lK + n~] D K~ . For the suppressed decay B~ — > [K + tt~]d^~ , both topology and 
phenomenology are similar to B~ — ► [K + tt~]dK~; our previous publication reported the 
first observation of this mode [7] . 




FIG. 1: Br -» [K + tt-] d K- and B -» [K + ir-} D ir- decays. 



ANALYSIS 

In this paper, we report an updated analysis of the suppressed decays B ± — > [K Z!? 'k ± ]dK ± 
and B^ — > [K^tt^dtv . In addition, the allowed decays 5 ± — ► [K ± tt t ]£ ) K ± and 
B ± — > [K ± -n T }£>'K ± are used as control samples to reduce systematic uncertainties. The 
same selection criteria for the suppressed decay modes are applied to the control samples 
whenever possible. The main changes with respect to our previous publication [7] are the in- 
clusion of additional data corresponding to 111 million BB pairs, and improved suppression 
of the dominant continuum background. Throughout this report, charge conjugate states 
are implied except where explicitly mentioned and we denote the analysed decay modes as 
follows: 

Suppressed decay B~ — > lK + n~] D h~ : B~ — > D snp h~ 
Allowed decay B~ -»■ [K-n + ] D h~ : Br -> D iaN hT (h = K, tt). 

The results are based on a data sample containing 386 million BB pairs collected at the 
T(45) resonance with the Belle detector at the KEKB asymmetric energy e + e~ collider [8]. 
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The Belle detector is a large-solid-angle magnetic spectrometer that consists of a silicon 
vertex detector (SVD), a 50-layer central drift chamber (CDC), an array of aerogel threshold 
Cerenkov counters (ACC), a barrel-like arrangement of time-of-flight scintillation counters 
(TOF), and an electromagnetic calorimeter (ECL) comprised of CsI(Tl) crystals located 
inside a superconducting solenoid coil that provides a 1.5 T magnetic field. An iron flux- 
return located outside of the coil is instrumented to detect K\ mesons and to identify muons 
(KLM). The detector is described in detail elsewhere [9]. Two different inner detector 
configurations were used. For the first sample of 152 million BB pairs, a 2.0 cm radius 
beampipe and a 3-layer silicon vertex detector were used; for the latter 234 million BB pairs, 
a 1.5 cm radius beampipe, a 4-layer silicon detector and a small-cell inner drift chamber were 
used [10]. 

Event selection 

D mesons are reconstructed by combining two oppositely charged tracks. These charged 
tracks are required to have a point of closest approach to the beam line within ±5 mm of 
the interaction point in the direction perpendicular to the beam axis (dr) and ±5 cm in 
the direction antiparallel to the positron beam axis (dz). A K/tt likelihood ratio P(K/n) = 
£>k/{£>k + £tt) is formed for each track, where Ck and C n are kaon and pion likelihoods, 
calculated using dE/dx measurements from the CDC, Cerenkov light yields in the ACC and 
timing information from the TOF. We used the particle identification requirement PiKj 7r) > 
0.4 and P(K/tt) < 0.7 for kaons and pions from D — > K7r decays, respectively. D candidates 
are required to have an invariant mass within ±2.5cx of the nominal D° mass: 1.850 GeV/c 2 < 
M(Ktv) < 1.879 GeV/c 2 . To improve the momentum determination, tracks from the D 
candidate are refitted according to the nominal D° mass hypothesis and the reconstructed 
vertex position (a mass-and- vertex-constrained fit). 

B mesons are reconstructed by combining D candidates with primary charged hadron 
candidates. For the primary charged tracks, we require P(K/tt) > 0.6 for the kaon in 
B~ — > DK~ and P(K/n) < 0.2 for the pion in B~ — > Dn~. The signal is identified 
by two kinematic variables, the energy difference AE = E D + E K -^-) — E heam and the 

beam-energy-constrained mass M bc = \J~E%~ a — (po + Pk-(tt-)) 2 , where E D is the energy 
of the D candidate, E^-u-) is the energy of the K~(n~) and i?beam is the beam energy, 
in the centre of mass (cm) frame. f>o and Pk-(-k-) ar e the momenta of the D and K~(ir~) 
in the cm frame. We define the signal region as 5.27 GeV/c 2 < M bc < 5.29 GeV/c 2 and 
—0.05 GeV < AE < 0.05 GeV. In the case of multiple candidates, which occurs in 1-2% 
of events with at least one candidate, we choose the best candidate on the basis of a \ 2 
determined from the difference between the measured and nominal values of M^ c . 

qq continuum suppression 

To suppress the large background from the two-jet like e + e~ — > qq (q = u, d, s, c) con- 
tinuum processes, variables that characterise the event topology are used. We construct a 
Fisher discriminant [11] of modified Fox- Wolfram moments [12], which we denote SFW. The 
Fisher coefficients are optimized by maximising the separation between signal events and 
continuum events. Furthermore, cos#s, the angle in the cm system of the B flight direction 
with respect to the beam axis is also used to distinguish BB events from continuum events. 
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These two independent variables, SFW and cos 6b, are combined to form a likelihood ratio 



^sig/ (^sig ^cont) 

r rSFW v pcos6 B 

*-sig(cont) — z -'sig(cont) X Z -'sig(cont) ' 

where £ sig and £ con t are likelihoods defined from SFW and cos9 B distributions for signal 
and continuum backgrounds, respectively. We optimize the TZ requirement by maximising 
S/y/S + N, where S and N denote the expected number of signal and background events in 
the signal region. The signal expectations are calculated from our previous results [7], the 
efficiency obtained from Monte Carlo simulation (given later), and the number of BB pairs; 
the background expectations are obtained using events in the M hc sideband (5.240 GeV/c 2 < 
Mb c < 5.265 GeV/c 2 ), with the extrapolation into the signal region based on Monte Carlo. 
For B- -> D sup K-(tt-) we require K > 0.90 (0.74), which retains 40.0% (65.7%) of the 
signal and removes 99.0% (94.3%) of the continuum background. 



Peaking backgrounds 

For B~ — > D sup K~ , one can have a contribution from B~ — > D°tt~ , D° — > K + K~ , which 
has the same final state and can peak under the signal. In order to reject these events, 
we veto events that satisfy 1.843 GeV/c 2 < M(KK) < 1.894 GeV/c 2 . The allowed decay 
B~ — > Df av /i~ can also cause a peaking background for the suppressed decay modes due to 
Kit misidentification. Therefore, we veto events for which the invariant mass of the Kn pair 
is inside the D mass cut window when the mass assignments are exchanged. Furthermore, 
three-body charmless decays B~ — * K + K~tt~ and B~ — > K + tt~tt~ can peak inside the signal 
region for B~ — > D sup K~ and B~ — > /} sup 7r~, respectively. These peaking backgrounds are 
estimated from the AE distributions of events in a D mass sideband, corresponding to 
±(2.5 - 10)<7 away from the nominal D mass (1.807 GeV/c 2 < M(Ktt) < 1.850 GeV/c 2 and 
1.879 GeV/c 2 < M(Kix) < 1.937 GeV/c 2 ). We fit these distributions, which are shown in 
Fig. 2, using a procedure similar to that used for candidate signal events (described later). 
For B~ — > D sup 7i~ , the peaking background estimated by fitting the plot is consistent with 
zero. Since the Standard Model prediction for the B~ — > K + n~n~ branching fraction is 
smaller than 10" 11 [13], this background contribution is ignored. On the other hand, for 
B~ — > D sup K~ , the peaking background yield in the D mass sideband is 7.3 t^o events, 
from which we expect 2.4i 2 ;o peaking background events inside the AE signal region. 

After applying all the cuts, the signal efficiencies are 14.6% and 24.5% for B~ — > D sup K~ 
and B~ — > D sup ir~, respectively. The signal yields are extracted by fitting the AE distribu- 
tions. 



Fitting the AE distributions 

Backgrounds from decays such as B~ — > Dp~ and B~ — > D*n~ are distributed in the 
negative AE region and make a small contribution to the signal region. The shape of 
this BB background is modelled as a smoothed histogram from generic Monte Carlo (MC) 
samples. The continuum background populates the entire AE region. We model its shape 
with as a first order polynomial. The signal AE distribution is modelled as the sum of two 
Gaussian distributions with a common mean. 
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FIG. 2: AE distributions for events in the D° mass sideband for (left) B — ► D sup K and (right) 
B~ — ► D sup ir~ . The signal shapes are modelled using the results of the B~ — > Df av h~ (h = K,ir) 
fit. 



TABLE I: Efficiencies and signal yields. For the B — > D sup K signal yield, the second value is 
after subtraction of the peaking background. 



Mode 


Efficiency (%) 


Signal Yield 


£- -► 


D slip K~ 


14.6 ±0.2 


2.4112/ o.o 


B- -► 


D sup ir~ 


24.5 ±0.3 


5011? 


5- -► 




14.5 ±0.2 


634l| 


B- -► 




24.9 ± 0.3 


14518 ± 125 



In the fit to the AE distribution of B~ — > D {av 7T~ , the free parameters are the position, 
widths, area and fraction in the tail of the signal peak, the slope and normalisation of the 
continuum component and the normalisation of the BB background. For the B~ — > Dfa v K~ 
fit, there is an additional component due to feed-across from D iav n~ , which we model with 
a Gaussian shape that has different widths on the left and right sides of the peak, since the 
shift caused by wrong mass assignment makes the shape asymmetric. The normalisation 
and shape parameters of this function are free parameters of the fit (all parameters which 
are floated in the B~ — > D lav it~ fit are again free parameters). 

For B~ — > D sup K~ and B~ — > D sup ir~, the signal and BB background shapes are mod- 
elled using the results of the fits to the corresponding favoured modes. The free parameters 
are the normalisations of the three components, and the slope of the continuum. The amount 
of feed-across from D snp n~ to D sup K~ is a free parameter, but is found to be negligible, as 
expected. The fit results are shown in Fig. 3. The numbers of events for B~ — > D sup h~ and 
D iav h~ are given in Table I. 
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FIG. 3: AE fit results for (top left) B~ -> D sup R-, (top right) -> -D S u P vr _ , (bottom left) 
— ► Df av K~, (bottom right) B~ — > Df av 7r~. Charge conjugate modes are included. 



RESULTS 

Ratio of branching fractions 

We calculate ratios of product branching fractions, defined as 



R 



Dh 



B{B~ 



B{B- 



D fav h 



Np sup h-/£D anp h- 



(h = K,n), 



where No snp h (^D fav /i) and eo sup h- ( e D !av h-) are the number of signal events and the recon- 
struction efficiency for the decay B~ — > D sup h~ (B~ — > D {ax h~), and are given in Table I. 
We obtain 

R DK = (O.Ol^(stat) ± l.O(syst)) x 1CT 3 , 
R Dn = (3.5 lg;f (stat) ±0.3(syst)) x 1(T 3 . 

Since the signal for B~ — > D sup K~ is not significant, we set an upper limit at the 90% 
confidence level (C.L.) of Rdk < 13.9 x 1CT 3 , where we take the likelihood function as a 
Gaussian distribution with width given by the quadratic sum of statistical and systematic 
errors, and the area is normalised in the physical region of positive branching fraction. 
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Most of the systematic uncertainties from the detection efficiencies and the particle iden- 
tification cancel when taking the ratios, since the kinematics of the B~ — > D sup h~ and 
B~ — > Df av h~ processes are similar. The systematic errors are due to uncertainties in the 
yield extraction and the efficiency difference between B~ — > D snp h~ and i?~ — > Df aY h~. The 
uncertainties in the signal shapes and the gg background shapes are determined by varying 
the shape of the fitting function by ±lcr. The uncertainties in the BB background shapes are 
determined by fitting the AE distribution in the region —0.07 GeV < AE < 0.20 GeV ignor- 
ing the BB background contributions — this is the largest source of uncertainty: the signal 
yields are affected by 7.4% and 28.4% for D sup 7i and D snp K (before peaking background sub- 
traction) respectively. The uncertainties in the efficiency differences are determined using 
signal MC. The total systematic error is the sum in quadrature of the above uncertainties. 

The ratio Rdk is related to 3 by 



R 



DK 



r l + r D + 2r B r D cos 3 cos 5, 



where [14] 



- 



A(B~ - 


-> D°K-) 




A(B~ - 


-> D°K~) 


1 


A(D° - 


-> K+TT-) 




A(D° - 


-> K-TT + ) 





5 = 5 B + 6d, 



0.060 ±0.003. 



and 8b ($d) is the strong phase difference between the two B (D) decay amplitudes. Using 
the above result, we obtain a limit onr B - The least restrictive limit is obtained allowing ±2cr 
variation on td and assuming maximal interference (03 = 0°, 5 = 180° or 03 = 180°, 5 = 0°) 
and is found to be r B < 0.18 at the 90% confidence level, as shown in Fig. 4. 




Belle Dalitz analysis 



FIG. 4: Constraint on tb from Rdk- 
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FIG. 5: AE fit results for (left) B -► D sup Tr~ and (right) B+ -► D sup ir+ . 



CP asymmetry 

We search for CP violating asymmetry in the 5 ± — > -D^pTr 1 * 1 mode. We fit the B + and 
B~ yields separately, and determine Ad-k as 

. B{B- -> AupTT-) - £(£?+ -> gsupTT+j 

D " - AupTT-) + B(B+ - L> sup 7T+) ' 

The fit results are shown in Fig. 5. We find 25.6 Ig^ — > -D S up7r~ events and 21.0 ly^ 
£> + — > £ ) sup 7r+ events, giving an asymmetry of 

A Dw = 0.10 ± 0.22 (stat) ±0.02 (syst), 

where systematic uncertainties arise from possible detector charge asymmetry (0.017; de- 
termined from the B^ — > D^tt^ sample [15]), and the B + and B~ yield extraction (0.016; 
determined as for Rd-k)- The total systematic error is obtained by taking the quadratic sum. 
The measured partial rate asymmetry Ad-k is consistent with zero. 

SUMMARY 

Using 386 million BB pairs collected with the Belle detector, we report studies of the 
suppressed decays B~ — > D sup h~ (h = K,ir). We do not observe a signal for B~ — > D snp K~, 
and place a limit on the ratio of B decay amplitudes re < 0.18 at the 90% confidence level. 
We have measured the CP asymmetry in the related mode, Ad-k = 0.10 ± 0.22 (stat) ± 
0.02 (syst). 
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